Field quantization by means of a single harmonic oscillator 



Marek Czachor 
Katedra Fizyki Teoretycznej i Metod Matematycznych 
Pohtechnika Gdanska, ul. Narutowicza 11/12, 80-952 Gdansk, Poland 

and 

Arnold Sommerferld Institiit fur Mathematische Physik 
Technische Universitat Clausthal, 38678 Clausthal-Zellerfeld, Germany 

A new scheme of field quantization is proposed. Instead of associating with different frequencies 
different oscillators we begin with a single oscillator that can exist in a superposition of different 
frequencies. The idea is applied to the electromagnetic radiation field. Using the standard Dirac- 
type mode-quantization of the electromagnetic field we obtain several standard properties such as 
^\ • coherent states or spontaneous and stimulated emission. As opposed to the standard approach the 

^\ ' vacuum energy is finite and does not have to be removed by any ad hoc procedure. 

0\ ' 

Q \ I. HARMONIC OSCILLATOR IN SUPERPOSITION OF FREQUENCIES 

<D 

. The standard quantization of a harmonic oscillator is based on quantization of p and q but w is a parameter. To 
^vq ' have, say, two different frequencies one has to consider two independent oscillators. On the other hand, it is evident 
that there exist oscillators which are in a superposition of different frequencies. The example is an oscillator wave 
packet associated with distribution of center-of-mass momenta. 

This simple observation raises the question of the role of superpositions of frequencies for a description of a single 
harmonic oscillator. We know that frequency is typically associated with an eigenvalue of some Hamiltonian or, which 
is basically the same, with boundary conditions. A natural way of incorporating different frequencies into a single 
harmonic oscillator is by means of the frequency operator 
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where all ujk > 0. For simplicity we have limited the discussion to the discrete spectrum but it is useful to include 
O ' from the outset the possibility of degeneracies. The corresponding Hamiltonian is defined by 
I 

Ch ; H = nn<g>- (a f a + aa^) (2) 
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^1 where a = X^^Lo V 71 + 1|ti)(ti + 1|- The eigenstates of H are \uj k , jk, n) and satisfy 

H\uj k Jk,n) = hw k (n+ ^j\^k, jk,n). (3) 

The standard case of the oscillator whose frequency is just ui coresponds either to O = ul or to the subspace spanned 
by \u> k , with fixed u>k = to. Introducing the operators 
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we find that 



The algebra of the oscillator is 
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[ a u> k ,j k ,a>l ltjl ] = ^ h u,i5j hjl \uj k ,j k )(w k ,j k \ <g> 1 (6) 

a ^k jk ji = ^k^i djkji ( a ^>k jk ) CO 

"l^Auji = 5 "k«>, 5 Mi( a tk,j k ) 2 - (8) 
The dynamics in the Schrodinger picture is given by 
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ihdt | V)=H\9) = Ml (g) (a f a+i (9) 

In the Heisenberg picture we obtain the important formula 

^ k , jk (t)=e iHt/H ^J h e- iHt/h (10) 

= \u k ,j k )(u k ,j k \ ® e-^'a = e-^a^O). (11) 

Taking a general state 

IV') = ^2 ^(^k,jk,n)\uj k , j k )\n) (12) 

we find that the average energy of the oscillator is 

(H) = {^m) = E \^kJk,n)\ 2 TuJ k (n+^j. (13) 

The average clearly looks as an average energy of an ensemble of different and independent oscillators. The ground 
state of the ensemble, i.e. the one with ip(u> k ,j k , n > 0) = has energy 

(H) = \ XJ \^k,J k ,0)\ 2 nLo k < oo. (14) 

The result is not surprising but still quite remarkable if one thinks of the problem of field quantization. 

The very idea of quantizing the electromagnetic field, as put forward by Born, Heisenberg, Jordan []l] and Dirac ||, 
is based on the observation that the mode decomposition of the electromagnetic energy is analogous to the energy of 
an ensemble of independent harmonic oscillators. In 1925, after the work of Heisenberg, it was clear what to do: One 
had to replace each classical oscillator by a quantum one. But since each oscillator had a definite frequency, to have 
an infinite number of different frequencies one needed an infinite number of oscillators. The price one payed for this 
assumption was the infinite energy of the electromagnetic vacuum. 

The infinity is regarded as an "easy" one since one can get rid of it by redefining the Hamiltonian and removing 
the infinite term. The result looks correct and many properties typical of a quantum harmonic oscillator are indeed 
observed in electromagnetic field. However, once we remove the infinite term by the procedure of "normal reordering" 
the resulting Hamiltonian is no longer physically equivalent to the one of the harmonic oscillators. For a single 
oscillator we can indeed add any finite number and the new Hamiltonian will describe the same physics. But having 
two or more such oscillators we cannot remove the ground state energies by a single shift of energy: Each oscillator has 
to be shifted by a different number and, accordingly, we change the energy differences between the levels of the global 
Hamiltonian describing the multi-oscillator system. And this is not just "shifting the origin of the energy scale" . 
Alternatively, one can add up all the ground state corrections and remove the overall energy shift by a different choice 
of the origin of the energy scale. This would have been acceptable if the shift were finite. Subtraction of infinite 
terms is in mathematics as forbidden as division by zero. (Example: l + oo = 2 + oo=>l = 2isas justified as 
l-0 = 2-0=*l = 2.) 

The oscillator which can exist in superpositions of different frequencies is a natural candidate as a starting point 
for Dirac-type field quantization. We do not need to remove the ground state energy since in the Hilbert space of 
physical states the correction is finite. The question we have to understand is whether one can obtain the well known 
quantum properties of the radiation field by this type of quantization. 

II. FIELD OPERATORS: FREE MAXWELL FIELDS 

The energy and momentum operators of the field are defined in analogy to H from the previous section 

H = Ku\\s,k\)(s,k\\ (g) ^(a^a + aa^ (15) 

s,k a 
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P = ^ ftK A |s,K A )(s,K A | <x> ^(a t a + aa t ) (17) 
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\ Yl M ^( a t,K X a s,K X + a s , Kx at tK3 ^ (18) 
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where s = ±1 corresponds to circular polarizations. Denote P = (H/c, P) and P ■ x = Ht — P ■ x. We employ the 
standard Dirac-type definitions for mode quantization in volume V 



A(t, £) = £ A/ — - (a^ x e-^e s ^- 3 + at^e^'e^e"^) (19) 
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= e iP^/^Ae' lP - x/n (20) 



k{t,x) = i£ J^fa^e-^e^e^ < KA e^ V^X% A ) (21) 



= e iP^I^Ee- lP - xln (22) 
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= e }P-*I^Be- ip - x l h . (25) 
Now take a state (say, in the Heisenberg picture) 

I*) = Yl *s,-A,«l s >«A,n) (26) 
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where |a s ,K A ) form a family of coherent states: 

a\a s ^ x ) = a s ^ x \a s ^ x ) (28) 

The averages of the field operators are 



(H>\A(t, = £ |*.,a A I 2 J (a SiKx e-^e s , Kx + a^e^e^) (29) 
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(*\E(t,x)\*) = \*sM x \ 2 J^(as^ x (0)e-^ x e s , Kx - < KA (0)e^'*e s * J (30) 
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These are just the classical fields. More precisely, the fields look like averages of monochromatic coherent states with 
probabilities |$ Sj £ A | 2 . The energy-momentum operators satisfy also the standard relations 



l - j <P~x(E{t,x)- E(t,x)+B(t,x)- B{t,xfj, (32) 

/ d 3 xE(t,x) x B(t,x). (33) 
Jv 
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It should be stressed, however, that these relations have a completely different mathematical origin than in the 
usual formalism where the integrals are necessary in order to make plane waves into an orthonormal basis. Here 
orthogonality follows from the presence of the projectors in the definition of a SiKx and the integration in itself is trivial 
since 
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E(t, x) ■ E(t, x) + B(t, x) ■ B(t, x)=E-E + B-B (34) 
E(t, x) x B(t, x)=Ex B. (35) 

Therefore the role of the integral is simply to produce the factor V which cancels with 1/V arising from the term 
1/y/V occuring in the mode decomposition of the fields. To end this section let us note that 

<*|ff|tt) = ]T ^ A |$ S , KA | 2 (K KA | 2 + \) (36) 
= ]T ^a|$ s , Ka | 2 (|^ Ka | 2 + \)- (37) 

S,K X 

The contribution from the vacuum fluctuations is nonzero but finite. 



III. SPONTANEOUS AND STIMULATED EMISSION 



The next test we have to perform is to check the examples that were responsible for the success of Dirac's quan- 
tization in atomic physics. It is clear that no differences arc expected to occur for single- mode problems such as the 
Jaynes-Cummings model. In what follows we will therefore concentrate on spontaneous and stimulated emission from 
two-level atoms. 

Beginning with the dipolc and rotating wave approximations we arrive at the Hamiltonian 

H = ^huj a 3 + ^J2 hLJx ( a l,K X a sMx + a s,Zx a l,z x ) + huJ od Y {g s ,R x a s M x (j + + ^aj ^cr_) (38) 
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where du — ) is the matrix element of the dipole moment evaluated between the excited and ground states, 

and g s .K X — i\J 2hu> x v '^ s <K>> ' ^ nc Hamiltonian represents a two- level atom located at x = 0. 
The Hamiltonian in the interaction picture has the well known form 

Hj = hu;odJ2 (9s,*^»°-» ^a s ,t x a + + gUs'^'^Us-)- ( 39 ) 
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Consider the initial state 

l*(°)>= *s',Z X >,rn\s',K\',m,+) 
s' ,k x / .rn 

= ^-'.^.ol*', Z'o, 0, +> + *s',<,nW, C n, +). (40) 

s' : k' q s' 

The states corresponding to n = play a role of a vacuum. As a consequence the vacuum is not represented here by a 
unique vector, but rather by a subspace of the Hilbert space of states. It is also clear that the energy of this vacuum 
may be nonzero since no normal ordering of observables is necessary. 

Using the first-order time-dependent perturbative expansion we arrive at 

!*(*)> = l*(o)> 

e -*(u>o-u>* )t _ I 

*-rr coo — u>\ » 
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E e -i(u> a -u> Xrl )t _ 2 
* s ,K„,nV / n+l5*s \s,K n ,n+ 1,-). (41) 
Ll>0 — LO\ 

One recognizes here the well known contributions from spontaneous and stimulated emissions. It should be stressed 
that although the final result looks familiar, the mathematical details behind the calculation are different from what 
we are accustomed to. For example, instead of 
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a\ uRi \s, it, m) ~ |si, Ki, 1; s, it, m), (42) 
which would hold in the standard formalism for Hi ^ it, we get simply 

=0, (43) 

a consequence of at - a[ s = 0. 
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